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The burgeoning field of Bose-Einstein conden- 
sation in dilute alkali and hydrogen gases has 
stimulated a great deal of research into the sta- 
tistical physics of weakly interacting quantum 
degenerate systems 1 ' 2 . The recent experiments 
offer the possibility for exploring fundamental 
properties of low temperature physics in a very 
controllable and accessible way. One current 
goal of experimenters in this field is to observe 
superfluid-like behavior in these trapped Bose 
gases, analogous to persistent currents in super- 
fluid liquid helium, which flow without observable 
viscosity, and electric currents in superconduc- 
tors, which flow without observable resistance. 
These "super" properties of Bose-condensed sys- 
tems occur because the macroscopic occupa- 
tion of a quantized mode provides a stabilizing 
mechanism that inhibits decay due to thermal 
relaxation 3 . Here we solve the time-dependent 
Gross-Pitaevskii equation of motion of the con- 
densate involving two hyperfine atomic states and 
show how to generate, with extremely high fi- 
delity, topological modes such as vortices that 
open the door to the study of superfluidity in 
these new systems. Our approach is inspired 
by recent experiments investigating a trapped 
condensate with two strongly coupled internal 
states 4,5 . We show how the interplay between the 
internal and motional dynamics can be utilized to 
prepare the condensate in a variety of interesting 
configurations. 

Since 1995, when Bose-Einstein condensation in a di- 
lute atomic gas was first observed 6-8 , experimenters have 
sought a method to create a vortex in this system. In a 
typical experiment, around one million atoms are trapped 
in a magnetic harmonic potential and cooled below the 
critical temperature so that condensation occurs into the 
lowest energy quantized mode. In the usual case, this 
ground state has no circulation. One proposed scheme 
for preparing the condensate in a vortex mode 9-17 is to 
distort the confining potential and mechanically rotate 
the trap during the cooling process. In this way, the 
lowest energy mode may be engineered to be circulat- 
ing about the axis of symmetry. Such an approach is 
in direct analogy with experiments on vortices in super- 
fluid helium — the asymmetry of the harmonic trap for 
the atomic gas plays the role of surface roughness of a 
rotating vessel. Although conceptually this method ap- 
pears promising for vortex generation in a trapped gas, 
so far technical difficulties have precluded its successful 
implementation. 



Instead of having the system condense into a vortex 
mode, an alternative approach is to allow the atoms to 
condense into the usual ground state and then dynam- 
ically generate the vortex from the non-rotating con- 
densate. Several theoretical proposals have been made 
along these lines which utilize the interaction between the 
atoms and a specific laser field consisting of a beam of 
photons with non-zero orbital angular momentum 18-20 . 

The method we present here makes use of both of the 
techniques mentioned — mechanical rotation and the cou- 
pling of internal states using an electromagnetic field. It 
is motivated by the recent observations demonstrating 
the cyclic twisting and untwisting of the order parame- 
ter of a gaseous condensate of rubidium atoms 4 . These 
recent experiments have suggested the possibility for har- 
nessing the interplay between the internal and motional 
dynamics in order to fundamentally alter the topological 
structure of the order parameter of the condensate 5 . 
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FIG. 1. Method for creating a vortex, (a) The two traps 
are rotated in the xy-plane about the z axis at a frequency 
oo r . (b) Transitions are simultaneously driven between the 
two internal states at the effective Rabi frequency Q e ff- A 
vortex mode possessing one unit of angular momentum can 
be prepared if uo r ~ ^ e ff- 

We treat the internal structure of the condensed atoms 
as a two- level system as illustrated in Fig. 1. The two 
states are confined in separate axially symmetric har- 
monic oscillator potentials with the same trap frequency 
uoq. As shown in Fig. 1(a), the trap centers are spatially 
offset by a distance ro and are rotated about the sym- 
metry axis at a frequency of uo r . Simultaneously an elec- 
tromagnetic field is applied that couples the two internal 
atomic hyperfine states causing the atoms to coherently 
cycle between levels as illustrated in Fig. 1(b). There 
are two parameters that characterize the coupling: the 
detuning and the power. The detuning 5 denotes the mis- 
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match of the frequency of the coupling electromagnetic 
field to the frequency difference between the two inter- 
nal atomic states. The power is characterized by the Rabi 
frequency £7, which is the rate at which population would 
oscillate between the two states if S were zero. When S 
is larger than so that the drive is off resonance, the 
population oscillations have small amp litude an d occur 
at the effective Rabi frequency ^ e ff = + £ 2 - 

The dynamical evolution of the condensate state \ijj(t)) 
for a gas of atoms in two different hyperfine states is 
governed by a nonlinear Schrodinger equation, known as 
the Gross-Pitaevskii equation, generalized to treat the 
internal state coupling 5 

d h 
ih—\i/;(t)) = [H ®l + l®-(n& x + 5a z ) 

+ (1) 

where {l,d- x ,a y ,a z } is the set of standard Pauli spin op- 
erators and for clarity we have used the tensor product 
(g) to explicitly separate the spatial and and internal op- 
erators. The free Hamiltonian Hq describes the motional 
dynamics of the atoms in a stationary harmonic trap cen- 
tered on the symmetry axis 

h 2 1 

Ho(r) = - — V 2 + -mu 2 r 2 + U n(r), (2) 
2m 2 

where m is the atomic mass. The term that makes the 
system nonlinear is the mean-field interaction energy, 
which depends on the local density of condensate atoms 
n(r). The coefficient U$ = 47rh 2 a/m is proportional to 
the scattering length of a binary collision a. The effect of 
displacing the trap centers and rotating them about the 
symmetry axis is described by 

Hi(r, t) = ft[f(r) cos(u r t) + g(r) sin(o; r t)], (3) 

where k is a coupling coefficient and / (r) and g(r) are 
functional prefactors. The explicit form of H\ determines 
the symmetry of the quantum state being prepared. To 
create a vortex state with one unit of angular momentum 
k = mo^ro, f(r) = x and g(r) = y in Cartesian coordi- 
nates, corresponding to the scheme depicted in Fig. 1(a). 
More general forms for / and g will be considered later 
with Eq. (3) then providing the definition of a generalized 
n and uo r . 

Before solving this problem in detail, we first present 
an intuitive picture of the underlying physics. The vortex 
state we wish to create has unit circulation correspond- 
ing to a macroscopically occupied single particle wave 
function with quantized angular momentum (L z )/h = 1. 
In order to see why our scheme couples a non-rotating 
condensate to this state, consider the frame co-rotating 
with the trap centers at angular frequency oj r so that H\ 
becomes time-independent. The free Hamiltonian in the 
co-rotating frame is given by 21 Ho — uo r L z . The energy of 
the vortex with one unit of angular momentum is there- 
fore shifted by Tiu r in the rotating frame relative to its 



value in the lab frame. When this energy shift compen- 
sates for both the energy mismatch TtS of the internal 
coupling field and the small chemical potential difference 
between the vortex and non-rotating condensate, reso- 
nant transfer of population may take place. 

For this picture to be valid, the following inequalities 
must be satisfied 

U0 r > CJ , S > U0 r W fi e ff, (4) 

which also imply « S. The first inequality splits sig- 
nificantly the energies of states of different angular mo- 
menta in the rotating frame. This point is crucial for ex- 
perimental implementation since it allows the vortex to 
be generated rapidly and is a feature absent from previ- 
ous dynamic coupling schemes 18-20 . The weak coupling 
limit given by the second inequality allows the resonance 
condition oj r « Q e ff to select energetically the desired 
state with high fidelity. 
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FIG. 2. Dynamical evolution to a vortex. These are results 
of a numerical integration of Eq. ([!]), with the condensate 
initially in the non-rotating ground state and in the internal 
state |1). The coupling drive is turned on at time t = 0, 
and is turned off at time t = t s by setting both Q and ro 
to zero. The top graph shows the fractional population of 
atoms in the |2) internal state. The small-amplitude rapid 
oscillations correspond to the cycling between internal levels 
due to the off-resonant coupling. The gradual rise of this line 
is due to coupling from the ground state to the vortex mode 
caused by the drive Hi in Eq. (|^). The bottom graph shows 
the angular momentum of the |2) state, in units of Planck's 
constant Ti. The rise and fall of this curve corresponds to 
a rapid cycling of the |2) atoms between the non- rotating 
condensate and the vortex. Once during each Rabi cycle, the 
angular momentum approaches unity and at that time the 
1 2) state wave function approaches a pure vortex mode. The 
inset shows the maximum amplitude of population transfer 
to the vortex as a function of the trap rotation frequency o; r , 
with A = Q e ff — Wr- The various parameters used in this 
calculation are: uo = 10 Hz, 6 = 200 Hz, uj r = 205.4 Hz, 
N = 8 x 10 5 atoms, m is the mass of the 87 Rb atom, a = 5.5 
nm 23 , and Q — 50 Hz and ro = 1.7/xm for t <t s . 



2 



Even though we require a large value for uj r , one may 
not use a simple equivalent time-averaged potential for 
the effect of rotating the trap centers without including 
the internal atomic state dynamics in the coarse-graining, 
which occur on the same time scale. 

In order to visualize the system dynamics, we present 
results from a numerical integration of the Gross- 
Pitaevskii equation 22 . Here we treat the system in two 
dimensions in the x?/-plane perpendicular to the symme- 
try axis. In Fig. 2 we plot the fractional population and 
the angular momentum per atom of the |2) state as a 
function of time. We initialize the system with all of 
the atoms in the |1) internal level and in the mean- field 
ground state of the trap. There are two striking features 
in Fig. 2. Remarkably, even though we are driving the in- 
ternal states far from the atomic resonance, a significant 
fraction of the population gets transferred to the |2) state 
over a time long compared to the fast Rabi oscillations at 
Cleft. Furthermore, we see that the angular momentum 
of state 1 2) oscillates rapidly at the frequency Q e ff. 




FIG. 3. Unit vortex preparation. Shown are the calculated 
densities and phases of the two states |1) and |2), at time 
t = 200 ms indicated in Fig. 2. At this time, about one third 
of the atoms are in the |2) state, which is in a pure vortex 
mode with unit angular momentum. A characteristic feature 
of a vortex mode is the 2nir phase wrap about the core, where 
n is an integer that is equal to unity in this case. An attractive 
feature of this preparation scheme is that the |1) atoms resid- 
ing in the core region provide a natural pinning mechanism 
for the vortex due to the mean-field repulsion. 

Our key idea is that by turning off the coupling at a 
precise time, t = t S: on a given Rabi cycle, the |2) state 
can be prepared to have unit angular momentum. The 
maximum possible population transfer to the vortex state 
using this scheme obeys a Lorentzian response curve as 
uj r is varied near exhibiting a narrow resonance. This 
is shown in the inset of Fig. 2, where A = — uj r , and 
will be discussed in more detail later. In Fig. 3 we show 
a snapshot of the densities and phases of the two com- 



ponents at time t = 200 ms. The snapshot illustrates the 
preparation of a high-quality vortex in state |2), with the 
|1) state providing a natural "pinning" mechanism that 
stabilizes the vortex core by providing a repulsive mean- 
field barrier along the symmetry axis. Consequently the 
size of the vortex core is determined by the spatial density 
profile of the non-rotating state, and not by the natural 
healing length. 

An intriguing property of our novel state preparation 
scheme is that the direction of circulation of the vor- 
tex can be opposite to that of the rotating trap centers; 
changing the sign of the detuning S with the direction 
of rotation of the trap centers fixed causes the vortex 
to rotate in the opposite direction. This is most easily 
seen by again considering the frame co-rotating with the 
potentials at frequency uo r . Vortices with opposite circu- 
lations experience opposite energy shifts in transforming 
to the rotating frame and therefore require opposite signs 
of detuning in order to achieve resonant coupling. 

Although the numerical calculations we present here 
are full solutions of Eq. (|l|) , an accurate and greatly sim- 
plified model may be derived using the separation of time 
scales given by the inequalities in Eq. (||) . This allows us 
to gain insight into the physical mechanisms giving rise 
to the behavior of this system. If we coarse-grain over 
rapidly oscillating terms in order to study the long time 
scale behavior, we arrive at an approximate solution for 
the system dynamics, which has the simple form 24 

\rl>(t)) = [a(t)co(t) |0 O ) + b(t)c n (t) |0 n )]|l) 

+ [b(t)c (t) |0 O ) + a*(t)c n (t) |0 n >]|2>. (5) 

Here |0o) and \(j) n ) are the spatial parts of the state, with 
|0o) being the ground state and \<j) n ) being the target 
state — the unit vortex for the case considered in Fig. 3, 
determined by the self-consistent lowest energy solutions 
of the two-component Gross-Pitaevskii equation 5 . For 
given populations of the two atomic states, the energy is 
minimized using a steepest descents algorithm with the 
constraint that |0o) is the nodeless ground state and \(j) n ) 
possesses the symmetry imposed by Hi, which for the 
unit vortex is a 2tt topological phase wrap on any closed 
loop containing the core. The rapidly varying coefficients 
in Eq. (||) are found to be given by a(t) = cos(fi e ff£/2) — 
i(5/Q e fi)sm(Q e tft/2) and b(t) — — i(Q/Q e fi) sm(Q e fit/2). 
The slow coupling of the spatial states is given by Co and 
c n , which satisfy 

■ h ( c N _ if -On - e - hA) f3{(j) n \f + ig\(/)o) \/c \ 
\c n ) ~ \ f3(fo\f ~ i9\M -eo-hA) J \ c n J ' 

(6) 

Here /3 = k,Q/\5\ is the coupling coefficient, / and g are 
defined in Eq. (||), and eo and e n are the energy eigenval- 
ues of |0o ) and |0 n ), respectively. The time scale for state 
preparation is determined by f3 and the overlap matrix 
element (cj) n \f + ig\(j)o). These two parameters also deter- 
mine the width of the response Lorentzian that is shown 
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in the inset of Fig. 2 with the maximum occurring not at 
zero frequency but displaced by an amount according to 
the splitting between the eigenenergies e n — eo; typically 
a small fraction of c^o- 

The symmetry of the coupling is determined by / + ig, 
which has the form x+iy as previously stated for the case 
of the unit vortex. In order to produce a vortex with n 
units of angular momentum (i.e. (L z )/Ti = n), / + ig 
should take on the form (x + iy) n . So for n = 2, we 
must construct an Hi in which / = x 2 — y 2 and g = 2xy, 
corresponding to a saddle potential rotating at uo r /2. In 
Fig. 4, we illustrate vortex generation with two and three 
units of angular momentum. 

U =311 




FIG. 4. Double and triple vortex preparation. Shown are 
the densities and phases modulo 2n of the \2) state after a dy- 
namical evolution similar to that shown in Fig. 2, but with dif- 
ferent symmetries of the drive Hi. In the case of (L z )/H = 2, 
f — x 2 — y 2 and g — 2xy, while for {L z )/h = 3, / = x 3 — 3xy 2 
and g = 3x 2 y — y 3 . In both cases, the system was evolved from 
the same initial condition as that for the calculation described 
in Fig. 2, with about one third of the atoms in the |2) state 
at the time t s . The values taken for the various parameters 
were the same, except for the trap rotation frequency, which 
was co r = 204.3 Hz for the (L z )/% = 2 case, and co r = 200.2 
Hz for the (L z )/h = 3 case. 

The generalization of our scheme for the preparation 
of macroscopic quantum states of arbitrary symmetry 
is straightforward. For example, to generate a mode 
with 



and higher angular momentum through collisional relax- 
ation. Even after the vortex state has been generated 
and the drive is turned off, collisions will be important 
in determining the stability of persistent currents where 
the topological structure of the two component order pa- 
rameter is known to be crucial 3 ' 12 . 

quadrupole 




FIG. 5. Dipole and quadrupole preparation. Different sym- 
metries were constructed using specific forms for the drive Hi 
in order to prepare the |2) state in non- circulating modes. In- 
stead of having a 2nn phase wrap corresponding to a current 
flowing around a central core, these modes have regions of 
constant phase separated by a discontinuous jump of ty where 
the wave function changes sign. To generate the dipole mode, 
we used f — x and g — with uo r — 205.4 Hz; while for the 
quadrupole, we used f — xy and g — with u r — 204.3 Hz. 

Note added: While this paper was in press, vortices 
in a two-component Bose-Einstein condensate have been 
created and observed using this scheme 25 . 
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righ t circulating vortex so that J — Re [x I iy) and 

g = 0. A quadrapole symmetry would be generated by 
/ = Im{(x + iy) 2 } and g = 0. The dynamical state 
preparation of these two examples is illustrated in Fig. 5. 

We do not anticipate collisions to play a significant 
role on the short time scale for vortex generation shown 
in Fig. 2 and we have implicitly neglected them in our 
zero temperature theory. However, the system is thermo- 
dynamically unstable on long time scales if the rotating 
drive Hi is left on, since atoms can access states of higher 
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